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We report on the first search for atmospheric and for diffuse astrophysical neutrino-induced show- 
ers (cascades) in the IceCube detector using 257 days of data collected in the year 2007-2008 with 22 
strings active. A total of 14 events with energies above 16TeV remained after event selections in the 
diffuse analysis, with an expected total background contribution of 8.3 ±3.6. At 90% confidence we 
set an upper limit of -E 2 $go%C£ < 3.6 x 10 -7 GeV ■ cm -2 • s _1 ■ sr -1 on the diffuse flux of neutrinos 
of all flavors in the energy range between 24TeV and 6.6 PeV assuming that "I> oc E~ 2 and that the 
flavor composition of the v e : i/ M : v T flux is 1 : 1 : 1 at the Earth. The atmospheric neutrino analysis 
was optimized for lower energies. A total of 12 events were observed with energies above 5 TeV. The 
observed number of events is consistent with the expected background, within the uncertainties. 



PACS numbers: 14.60.Lm, 95.85.Ry, 95.55.Vj 
I. INTRODUCTION 

The origin of high-energy cosmic-rays is an area of ac- 
tive research in astrophysics. The sites where cosmic rays 
are accelerated are expected to produce high energy neu- 
trinos. Many types of objects, ranging from supernovae 
and gamma-ray bursters to active galactic nuclei [1] , have 
been proposed as point sources of high energy neutrinos 
and many searches for such sources have been made [2J, 
yielding results consistent with background only assump- 
tions. If there are many point sources, each with an un- 
observably low flux, then the aggregate flux may still be 
observable as a diffuse flux. 

Diffuse searches rely on the energy spectrum of the de- 
tected events to separate an extraterrestrial signal from 
atmospheric neutrinos produced in the interaction of cos- 
mic rays with atomic nuclei in the Earth's atmosphere. 
Low energy (below ~ 10 GcV) atmospheric muon and 
electron neutrinos have been observed in underground de- 
tectors [3]. At higher energies, from 100 GeV to 400 TeV, 



neutrino telescopes have measured the spectrum of at- 
mospheric [4]. In this energy range, the flux of v e is 
expected to be lower by about a factor of 20 [5] and has 
not been observed. 

The main component of the atmospheric neutrino spec- 
trum is produced by the decays of tt and K ± . Asymp- 
totically it can be parametrized by dN v /dE l/ oc E~ 3 - 7 , 
where E v is the neutrino energy [6]. Decays of hadrons 
containing charm and bottom quarks form an additional 
component that is expected to be close to the primary 
cosmic-ray spectrum, dN v jdE v oc E~ 2 - 7 [7HS], and pro- 
duces nearly equal numbers of and v e . These prompt 
neutrinos are expected to dominate the v e spectrum at 
energies above ~ 30 TeV [S]. The production of v T is 
expected to be negligible. 

Fermi acceleration of charged particles in magnetic 
shocks followed by collisions with matter or radiation 
between the source and the Earth naturally leads to 
an energy spectrum for extraterrestrial neutrinos that 
is harder than that for atmospheric neutrinos, typically 
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close to dN v /dE v oc E v 2 . This allows diffuse extrater- 
restrial neutrinos to be visible as a hard component to 



the observed spectrum. The ratio of the v e 



v T flux 



in a single astrophysical source depends on the neutrino 
energy [10] - At moderate (high) energies, the neutrino 
flavor flux ratio behaves like the one from a pure pion 
(muon-damped) source, leading to an observed 1:1:1 
(1:1.8:1.8) ratio at the Earth after taking into account 
neutrino oscillations. The energy at which a flavor ratio 
transition occurs thus depends on the properties of the 
source 10 . The neutrino flux is not known, although it 
is expected to be below the Waxman-Bahcall bound [TT] . 

Previous searches for a diffuse flux have been per- 
formed with muon neutrinos [6], and with cascades [12] 
[T3"] . Cascades are the particle showers (electromagnetic 
and hadronic) initiated by charged current interactions of 
v e and v T and the neutral-current neutrino interactions 
of all three flavors in a medium. In the charged-current 
interactions, an average of 80% of the (high energy) neu- 
trino energy goes into the produced lepton [14] . For is e , 
this leads to an electromagnetic shower, while for u T the 
character of the lepton-induced shower depends on the r 
decay mode. The remainder of the energy is transferred 
to the target nucleon, producing a hadronic cascade. In 
the neutral-current interactions, the neutrino transfers 
a fraction of its energy to the target nucleon producing 
only a hadronic cascade. A typical cascade deposits its 
electromagnetic energy in a thin cylinder about 30 cm 
in radius and 5 m in length. Hadronic energy is de- 
posited over a larger volume, about 11 m long and 75 
cm in diameter. IceCube observes the Cherenkov radi- 
ation produced by the charged secondary particles from 
neutrino-nucleon interactions through an optical sensor 
array. While a neutrino-induced muon has a track-like 
signature in IceCube, a cascade event looks effectively 
like a point source of Cherenkov light in the detector. 

For diffuse searches, cascades from all flavor v inter- 
actions have two advantages over tracks from in- 
teractions, despite their inherently poor angular reso- 
lution compared to muon tracks. The first is that the 
background from atmospheric neutrinos is lower than 
for v^. Second, because of their short shower length, 
the cascades are well-contained in the detector, with a 
Cherenkov light output proportional to the shower en- 
ergy, so the shower energy is well measured. The detec- 
tor acts as a calorimeter. Since the energy spectrum of 
extraterrestrial neutrinos is expected to be harder than 
the atmospheric neutrino spectrum, searching for a break 
in the energy spectrum with cascades is easier than with 
muons, both due to the expected break being at a lower 
energy in the cascade channel than the muon channel (a 
consequence of lower fluxes of atmospheric v e than v^), 
and better intrinsic energy resolution of cascades over 
muons. 

This paper reports on searches for diffuse extraterres- 
trial and for atmospheric neutrino-induced cascades us- 
ing 257 days (livetime) of data collected in the year 2007- 
2008 with a partially completed IceCube detector consist- 
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FIG. 1: (Color online) a) The filled circles show the posi- 
tions of the strings in the x — y (horizontal) plane for the 
22-string detector configuration. The lines show the bound- 
aries of the fiducial volume and are described in the text, b) 
The reconstructed center-of-gravity position x (COGX) ver- 
sus y (COGY) for events passing the cascade online filter. 
The right axis shows the rate [Hz]. 



ing of 22 of the planned 86 strings. The IceCube detector 
and data sample are described in section [TTJ. Section |III| 
describes the analysis. Results are given in section |IV| 
and a summary follows in section [Vj 



II. THE ICECUBE DETECTOR AND DATA 
SAMPLE 

The IceCube detector is composed of vertical strings 
of optical sensors which are deployed in the glacial ice 
at the South Pole. The sensors detect Cherenkov radi- 
ation from charged particles produced in neutrino inter- 
actions. The strings are deployed on a 125 m triangu- 
lar grid. Each string contains 60 digital optical modules 
(DOMs), mounted between 1450 and 2450 m below the 
surface (17 m spacing). Each DOM contains a 10-inch 
photomultiplier tube (PMT) [15] , and a data acquisition 
system in a pressure vessel. The DOMs are most sensi- 
tive to photons with wavelengths 300 to 600 nm. The 
Hamamatsu R7081-02 PMTs have a peak quantum effi- 
ciency of ~ 25% and are operated at a typical gain of 
10 7 . The PMT output is linear within 10% up to a cur- 
rent of 50 mA; at a gain of 10 7 , this corresponds to 31 
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detected photoelectrons per nanosecond. The DOM and 
string performance is described elsewhere \16\ 117] . 

The data acquisition system records the arrival times of 
the detected photoelectrons. It uses two waveform digiti- 
zation systems which record the arrival time with a time 
accuracy of about 2 ns and a wide dynamic range |16lll8j . 
One system records data at 300 megasamples/second 
(MSPS) for 400 ns after the first photon is detected with 
14 bits of dynamic range. The second only has 10 bits of 
dynamic range but records data for 6.4 /zs at 40 MSPS. 
The dead-time of the system is less than 1%. Each DOM 
is activated (launched) when a single photoelectron is 
detected, and the data are sent to the surface when two 
adjacent (nearest or next-to-nearest) neighbors record a 
hit within 1 fis. 

The data were collected between May 2007 and April 
2008 when IceCube consisted of 22 active strings with 
1320 DOMs. The detector configuration is shown 
schematically in Figure The point (0,0,0) is the cen- 
ter of the complete 86-string IceCube detector. The main 
physics trigger was a "simple multiplicity trigger" that 
required photon signals from at least 8 DOMs, with the 
additional requirement of accompanying hits in any of 
the two neighboring DOMs on a string, each above a 
threshold of 1/6 single photoelectron and within a 5 fxs 
coincidence window. The average trigger rate was ~ 550 
Hz, driven by atmospheric muons, and exhibited about 
±10% seasonal variation. Online filters were employed 
to preselect a data sample for satellite transmission, that 
was used for analysis. 

The online filter relevant to this paper used two first- 
guess reconstruction algorithms [19]. One algorithm as- 
sumed that all hits can be projected onto a line consistent 
with a particle traveling at a specific velocity. The second 
algorithm quantified the sphericity of the hit topology. It 
used the center-of-gravity (COG), defined as the photon 
signal amplitude weighted mean of all hit DOM positions, 
as a first-guess vertex position. The response of the algo- 
rithms was studied with Monte Carlo simulation. Online 
selection criteria were developed using these simulations 
to reject 97.5% of the background events, while retaining 
70% of the cascade signal events. The average rate after 
this filtering was 19 Hz. 



III. ANALYSIS 

Even after online filtering, the data are dominated by 
atmospheric muons produced in interactions of cosmic 
rays with nuclei in the Earth's atmosphere. It is thus nec- 
essary to develop stringent offline selections to reject this 
background while retaining cascade signal events. Two 
such analyses were performed independently. One anal- 
ysis focused on cascade events induced by high energy 
extraterrestrial neutrino interactions, while the other 
searched for lower energy atmospheric neutrino-induced 
cascade signals. The development of the selections for 
both analyses relies on using 10% of the recorded data 



while keeping the remaining 90% "blinded". When the 
selections have been fixed, the 10% sample is discarded 
and the remaining data is unblinded. The physics results 
were evaluated from the remaining 90% of the recorded 
data only after the analysis selections were finalized and 
were thus free of statistical bias. In addition, exten- 
sive Monte Carlo simulations of background and signal 
events, as well as the IceCube detector response, were 
used. These simulations used importance sampling and 
weighting techniques to overcome computing limitations 
and are described in section [III A| This is followed by a 
description of the two analyses in section 
atic uncertainties are described in section 



IIIB 



IIIC 



System- 



A. Monte Carlo Simulations 

The Monte Carlo generator ANIS [20] was used to gen- 
erate neutrinos of all flavors in the energy range from 
~ 10 GeV to 1 EeV at the surface of the Earth and to 
simulate subsequent neutral and charged current interac- 
tions in the Earth. The neutrino spectra were generated 
with energy distributions of E~ 2 and E~ l , and were re- 
weighted to conventional and prompt atmospheric neu- 
trino flux predictions or to an astrophysical neutrino flux 
of E 2 & = 1 x 10~ 6 GeV • cm~ 2 • s _1 • sr _1 , depending on 
the analysis goal. The CTEQ5 parton distribution func- 
tions |21j were used in evaluating the neutrino cross sec- 
tions. The Preliminary Reference Earth Model [22] was 
used as the density profile to simulate high energy neu- 
trino absorption. The simulations included the Glashow 
resonance v e + — > W~ — > X for D e energies near 6.3 
PeV [23]. 

CORSIKA [24] was used to simulate air showers from 
cosmic ray collisions in the Earth's atmosphere and ob- 
tain a sample of background atmospheric muons. The 
cosmic ray spectrum was based on the Horandel polygo- 
nato model [25 , and primaries up to iron were simulated 
in the energy range from 600 GeV up to 10 11 GeV. The 
SIBYLL parameterization [53] was used in CORSIKA to 
model hadronic interactions. The effects of the Earth's 
magnetic field were included in the simulations. By over- 
laying independent CORSIKA showers, we simulated the 
coincident muon background, which comes from multi- 
ple cosmic ray showers triggering the detector within its 
readout window. 

High energy background muons are likely to incur large 
radiative energy losses, which may mimic cascade sig- 
nals. To obtain an adequate sample of such high energy 
background muons in the simulations, it was necessary to 
impose additional selections in the CORSIKA event gen- 
erations. We have imposed threshold selection criteria on 
the primary cosmic ray energy and on the energy of indi- 
vidual muons, as was done in the AMANDA analysis [12] , 
to increase the effective livetime and Monte Carlo statis- 
tics for cosmic rays in the energy range from O(10 2 ) TeV 
to 0(1) PeV. The effective livetimes of these simulation 
samples are given, together with the threshold values, 
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TABLE I: The effective livetimes of the CORSfKA Monte 
Carlo simulation samples, generated with standard energy 
thresholds of 600 GeV for primary cosmic rays and 273 GeV 
for secondary muons and with high-energy thresholds of 
40 TeV and 5 TeV, respectively. 



Primary Cosmic Ray Energy 


Monte Carlo sample 




standard 


high-energy 




[days] 


[days] 


10 TeV 


1 




100 TeV 


10 


50 


1 PeV 


100 


500 


10 PeV 


1000 


5000 



in Table [T] for representative cosmic ray energies. The 
high-energy sample was used only to assess background. 

A Muon Monte Carlo simulator [STJ was used to propa- 
gate secondary muons through the ice. It simulates muon 
stochastic energy losses from ionization, bremsstrahlung, 
photo-nuclear interactions, and pair production. 

A Cascade Monte Carlo [28] was used to simulate 
the longitudinal development of electromagnetic and 
hadronic cascades. It also accounts for the lower 
Cherenkov light output from hadronic cascades compared 
to that from electromagnetic cascades. 

Cherenkov light emission and subsequent photon prop- 
agation through ice was simulated with the Photonics 
simulation package 29 . In these simulations, the opti- 
cal properties of the ice were described by a "calibrated 
ice model" . This model was constructed from extensive 
AMANDA measurements of light propagation in South 
Pole ice made with artificial in-situ light sources (pulsed 
and steady LED sources and nitrogen lasers) [30] . These 
measurements were largely decoupled from light source 
and detector characteristics by using timing information 
of detected single photoelectrons. They determined the 
relevant wavelength and depth dependences of the optical 
scattering and absorption lengths, down to the deepest 
AMANDA depths. Using the fact that scattering and 
absorption are highly correlated with the concentration 
of insoluble dust particles in the ice, the model was ex- 
trapolated to the greater IceCube depths (from 2100 m 
down to bedrock at about 2800 m) with dust concentra- 
tion data from an Antarctic ice core [31) . which were 
scaled to fit the measured scattering and absorption pa- 
rameters at AMANDA depths. The ice properties are 
thus less understood at depths greater than 2100 m. 

PMT response simulators were used for each DOM in 
the detector to relate light input and current output. 
The simulated currents were then propagated through 
response simulators of the digitization electronics, local 
coincidence signaling, and event triggering. The same 
processing and filtering was applied to simulated and 
recorded data. 



B. Event selections 

After online filtering and transfer, the data were passed 
through several maximum-likelihood based reconstruc- 
tion algorithms [THl 132] in order to suppress muon back- 
ground while retaining signal events. The algorithms as- 
sumed a single track, two track, and a point-like cascade 
hypothesis and obtained the corresponding vertex posi- 
tions and times for each event. Following the reconstruc- 
tion, a sequence of selection criteria was applied to the 
data. The retained data are identified by different levels, 
starting from Level-3 (Level-1 refers to triggered online 
data and Level-2 to the data before offline reconstruc- 
tion). 

At Level-3 only those events were retained for which 
the single muon track reconstruction found a zenith an- 
gle greater than 73°. In addition, a requirement was 
imposed on the difference of the log-likelihood, C, of 
the cascade and track reconstructions for each event, 
£(cascade) — £(track) > —16.2, to preferentially select 
cascade-like events. The selection criteria were chosen so 
that the signal efficiency was about 80% in both the at- 
mospheric and extraterrestrial neutrino searches and the 
muon background rate was reduced by a factor of ~ 5, 
with respect to the online filter. 

Both analyses used the same cascade energy recon- 
struction algorithm [33 which was applied to events pass- 
ing Level-3. This analytical cascade-energy calculation 
takes into account the variation of ice optical properties 
with depth: 

£rcco/GeV = ^ r. (1) 

The n pe is the observed number of photoelectrons in a 
DOM and [1q (jv/dom) is the number of photoelectrons 
expected at a given DOM position rboM from a 1 GeV 
cascade with a vertex position r v . The /io is taken from 
a table generated by the Photonics Monte Carlo pack- 
age [23]. In Eq. [I] the noise term has been neglected. 
It biases the reconstructed energy towards lower values 
and thus worsens the performance of the cascade energy 
reconstruction. 

The subsequent selections, described below, were opti- 
mized separately for the two analyses. 

Extraterrestrial event selections 

Energetic bremsstrahlung from muon tracks outside 
the detector and muon tracks that intersect only part 
of the detector can mimic an uncontained cascade sig- 
nal. To reject this type of background, events were re- 
quired to have a topology consistent with a cascade sig- 
nal that originated inside the IceCube instrumented vol- 
ume at Level-4. Specifically, the four earliest hits in the 
event were required to be inside the fiducial volume of 
the detector in the horizontal x-y coordinates, as de- 
picted in Figure [Ik (dotted lines). Events were also 
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rejected if any of the four earliest hits occurred in the 
eight topmost DOMs. Approximately 5% of the data 
and 7% of the muon background events from the COR- 
SIKA Monte Carlo remained with respect to the previous 
selection level. The difference between absolute rates in 
the data and the CORSIKA Monte Carlo is addressed in 
section IIII CI ANIS Monte Carlo simulations show that 
approximately 13% of the signal was retained at Level-4. 

At Levcl-5, events were required to have hits in at least 
20 DOMs. To reject events with multiple muon tracks 
within the IceCube read-out window, an upper limit of 
5/xs was set on the event duration, defined as the time 
difference between the last and the earliest hit in the 
event. A two-track reconstruction was performed on the 
retained events and the tracks were required to coincide 
to within 1 fis. The center-of-gravity of the hits and the 
vertex coordinates from the cascade reconstruction algo- 
rithm were required to coincide in x and y to within 60 m. 
Monte Carlo simulation studies indicate that these selec- 
tions reduce the coincident muon rate to less than 20% 
of the 8mHz expected background muon rate at level 5. 
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FIG. 3: (top) Normalized distribution of the reconstructed 
cascade energy after level-6 event selections for the data 
(filled circles), muon background Monte Carlo (continuous 
histogram) and signal neutrino (dashed histogram), (bottom) 
Absolute event rate versus primary Cosmic Ray energy after 
level-6 event selections with an additional selection criterium 
on the reconstructed cascade energy of 6.3 TeV for the stan- 
dard Corsika Monte Carlo (continuous histogram) and high- 
energy optimized Corsika Monte Carlo (open triangles). 
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FIG. 2: Normalized distributions of (top) zenith angle and 
(bottom) reduced log likelihood ratios after Level-5 event se- 
lections for the data (filled circles), muon background Monte 
Carlo (continuous histogram) and signal neutrino (dashed his- 
togram) . 

The smaller event sample allowed for the use of more 
CPU-intensive reconstruction algorithms. At Level-6, 
the single-track reconstruction was iterated 32 times with 
randomly selected seed tracks to ensure that the final re- 
constructed tracks did not originate from local minima. 



The zenith angle distribution is shown in the top part 
of Fig. [2] The selection on zenith angle was refined to 
reject events with a reevaluated single-track zenith angle 
smaller than 69°. The comparison of the event recon- 
struction probabilities with a single-track and a cascade 
hypothesis was revisited as well. Only those events for 
which the ratio of the reduced log-likelihood, r£, for the 
event reconstruction under a single-track and under a 
cascade hypothesis was greater than 0.95 were retained 
for further analysis. The r£ is defined as the negative 
log-likelihood normalized by the number of degrees of 
freedom. The normalized distribution of the track and 
cascade reduced likelihood ratio is shown in the bottom 
part of Fig. [2J The Monte Carlo describes the shapes 
of experimental data distributions very well. The rates 
observed in the data and obtained from simulations are 
given in Table [TT] and discussed in section |HI C| 

Figures [3] (top) and [4] show comparisons between 
the Monte Carlo and the data for the reconstructed en- 
ergy, and the cascade reconstructed vertex positions af- 
ter Level-6 selections. The shapes agree very well, ex- 
cept vertex position z at the largest depths of the de- 
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FIG. 4: Normalized distributions of cascade reconstructed vertex position: x component (left panel), y component (middle 
panel) and z component (right panel) after level-6 event selections for the data (filled circles), standard muon background 
Monte Carlo (continuous histogram). 



tector where the ice properties are not well described by 
our simulations. Figure [3] (bottom) shows the absolute 
event rate as a function of primary Cosmic Ray energy for 
the standard and high-energy optimized Corsika Monte 
Carlo after level-6 event selections with an additional se- 
lection criterium on the reconstructed cascade energy of 
6.3 TeV. Since the two spectra are consistent within un- 
certainties, at Level-7 we use the high-energy optimized 
Corsika Monte Carlo. 

In order to further distinguish between signal-like and 
muon-like events, at Level-7 we used the multiplicity and 
spatial distribution of the hit DOMs. For every event we 
calculated the distances between each hit DOM and the 
reconstructed cascade vertex. Event-by-event, the mean 
distance, D, was evaluated and used as the half-radius 
in determining fill-ratio F, defined as a fraction of hit 
DOMs within a sphere centered on the cascade vertex for 
the event. Monte Carlo studies of signal and background 
exhibit different correlations of D and F, as shown in Fig- 
ure [5j For signal-like events, the average F value (sphere 
density) increases as a function of D, while the opposite 
is true for the single muon events. Most of coincident 
muon events have very small values of F, independent of 
D. Figure [5] also shows the distribution of D and F in 
the data and Monte Carlo simulations. The background 
Monte Carlo is in good agreement with the data for the 
events with the cascade reconstructed vertex position 
z > —300m as shown in Figure [5] (bottom-right panel). 
This is not the case for the events with the cascade recon- 
structed vertex position z < —300 m, where the ice is less 
well understood and makes some muons look like non- 
contained cascades (spherical shape and high DOM mul- 
tiplicity). The muon background events were removed 
by increasing a threshold on the minimum number of hit 
DOMs to 60 and imposing a threshold E rcco > 16 TeV on 
the reconstructed cascade energy. These selections also 
suppress contributions from atmospheric neutrinos from 
pion and kaon decays. After all selections, four events in 
the muon background simulation samples and no events 
in the 10% of the data remained [53] . 




D[m] D[m] 



FIG. 5: (Color online) The fill-ratio, F, versus the distance, 
D, as defined in the text for (top-left panel) the Monte Carlo 
simulated signal sample, (top-right panel) the Monte Carlo 
simulated muon background sample, and the data at the bot- 
tom part of the detector (bottom-left panel) and at the top 
part of the detector (bottom-right panel). The right axis 
shows the rate [Hz]. The dashed lines show the selection 
boundaries applied at Level-7 in the analysis. 



Atmospheric Neutrino Event Selection 

The atmospheric neutrino analysis used artificial neu- 
ral networks implemented in the ROOT TMVA pack- 
age [35] to reject background. The neural networks had 
two hidden layers, the first with N + 1 neurons and the 
second with iV neurons, where N is the number of dis- 
criminating variables. 

At Level-4, five input variables were used; the track 
zenith angle of the 32-fold iterative muon reconstruction, 



TABLE II: Event rates at different selection levels for the extraterrestrial and atmospheric analyses, described in 
the text, for the data and for Monte Carlo simulations of atmospheric background muons and of atmospheric and 
astrophysical neutrinos. An astrophysical neutrino flux of -E 2< I > m odci = 1-0 x 10 _6 GeV • cm -2 • s _1 ■ sr _1 was used 
for all v flavors. The statistical uncertainties in the event rates are smaller than 10%. 
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the reduced likelihood parameter from the cascade ver- 
tex reconstruction, the number of DOMs that register 
a single photoelectron divided by the total number of 
photoelectrons seen by all DOMs, the number of direct, 
unscattered photons assuming a point source at the re- 
constructed cascade vertex, and the difference in z vertex 
positions when the cascade vertex is reconstructed using 
only the earliest registered hits and using only the latest 
registered hits for the so-called split cascade reconstruc- 
tions. Figure [6] shows the classification score, Qa, of 
the neural network trained from these variables. Only 
events with Qa > 0.4 were retained. In addition, the 
reconstructed cascade was required to have an energy 
larger than 2TeV and be contained within the detec- 
tor by imposing an upper value of 1.4 on the parallelo- 
gram distance, a, of its horizontal vertex position. The 
six innermost strings span a parallelogram in x and y, 
whose edges define a = 1. The neighboring strings span 
a stretched parallelogram with a = 2. Parallelograms 
of a = 1 (a = 2) are shown as dashed (dotted) lines in 
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FIG. 6: Normalized distributions of the Level-4 neural net- 
work classification score. 



Fig. [Tp). The containment selection criterium a < 1.4 
thus requires the vertex to lie at least 75 m from the edge 
of the detector. With these selections, background was 
reduced by a factor of 1700 and 2.4% of the signal events 
were retained. 

At Level-5, we chose to train two individual neural net- 
works separately because of limited remaining statistics 
in the Monte Carlo training samples and to achieve ade- 
quate performance. The first neural network was trained 
with the number of direct hits from the reconstructed 
cascade vertex, the total number of observed photoelec- 
trons in all DOMs, the difference between log-likelihoods 
from the cascade vertex reconstruction and from the 32- 
fold iterative muon track reconstruction, and the recon- 
structed cascade energy. These four variables correlate in 
a complex way that merits the use of a neural network. 
The second neural network was trained with six variables: 
the distance from the cascade vertex to the first hit in the 
event, the fill-ratio F, the reduced log-likelihood parame- 
ter from the cascade vertex reconstruction, the difference 
in z vertex positions for the two split cascade reconstruc- 
tions, the track zenith angle of the 32-fold iterative muon 
reconstruction, and a defined and used before at Level-4. 
Figures [7] and [8] show the distributions of selected input 
variables, whereas Fig. [9] shows the two neural network 
classification scores Qb and Qo 

The three neural network classification scores, Qa, Qb 
and Qc were optimized for E rcco > 5 TeV, and their 
product, Q* = Qa x Qb x Qc, was used at Level-5. In 
addition, events with a reconstructed cascade vertex in 
the topmost 60 meters of the detector were rejected. 
Multi-muon background from coincident cosmic ray air 
showers was efficiently rejected by a relatively loose selec- 
tion criterium on the reduced likelihood parameter from 
the cascade vertex reconstruction. The final selection was 
Q* > 0.73 for E ICCO > 5 TeV. 
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FIG. 7: Normalized distributions (after Level-4 atmospheric selections) of (left panel) reconstructed energy (middle panel) 
number of direct hits (right panel) total number of hits. 
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FIG. 8: Normalized distributions (after Level-4 atmospheric selections) of (left panel) distance from the cascade vertex to 
the first hit in the event (middle panel) difference in z vertex positions for the two split cascade reconstructions (right panel) 
fill-ratio from the mean hit distance. 
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FIG. 9: Normalized distributions of the neural network 
classification scores (top) Qb and (bottom) Qc. 



C. Normalization and Systematic Uncertainties 

The extraterrestrial and atmospheric analyses have 
common sources of systematic uncertainty. The largest 
contributions to our systematic uncertainty estimate 
arise from our limited knowledge of the optical proper- 
ties of the ice and from uncertainties in the cosmic-ray 
flux and composition. Other significant contributions re- 
sult from uncertainty in the DOM detection efficiency 
and from uncertainties in the neutrino cross-sections and 
the light output from the cascades. These sources affect 
signal and background estimates differently. We describe 
signal and background separately below. 



Muon Background Systematics 

CORSIKA Monte Carlo simulations adequately repro- 
duce the shape of the observed spectra in both analyses 
but systematically underestimate the absolute rate. In 
the extraterrestrial analysis, an empirical factor to nor- 
malize the rate of simulated background to the corre- 
sponding data was found to be 3 ± 1, where the uncer- 
tainty covers the difference between the data and COR- 
SIKA rates at selection levels up to 6 (see Table |rij ). 
Background events outnumber signal events by at least 
3 orders of magnitude in these event samples. In the 
atmospheric analysis, the same factor 3 ± 1 was found 
for Q* < 0.4, where any cascade signal is expected to 
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FIG. 10: The effective area, Ajff, for (left panel) u e , (middle panel) v^, and (right panel) v r versus neutrino energy at trigger level 
(dashed curves) and at analysis level after all selections for both the extraterretrial and the atmospheric analysis (continuous 
curves). 



be vanishingly small. The factor was the same, within 
uncertainties, at larger Q* values. 

Simulations based on an alternate ice model |36] were 
performed to assess the possible effects on the empirical 
normalization factor caused by incomplete description of 
the optical ice properties and necessary approximations 
in their implementation in the simulations. The alter- 
nate model is based on a global fit to recorded charges 
in very bright events generated by IceCube LED sources 
sampling all detector depths. The background rates sim- 
ulated with the alternate ice model differed from the rates 
observed in the data by a factor of two. This is consis- 
tent, to within the assigned uncertainties, with the dif- 
ference between the rates simulated using the calibrated 
ice model and those observed in the data. 

The remaining difference in simulated and observed 
rate is ascribed to the combined effect of uncertainties in 
the cosmic ray flux for protons in the energy range near 
the knee and in the absolute energy scale. 

Neutrino Systematics 

Simulations show that neutrino events are much less 
affected by uncertainty in the ice model. We estimate a 
20% uncertainty in the rates by comparing simulations 
with the calibrated and alternate ice models for both 
atmospheric and E~ 2 energy spectra. 

The uncertainty in the DOM sensitivity is taken as 8%, 
based on the measured uncertainty in the PMT sensitiv- 
ity |15j . For an E~ 2 neutrino energy spectrum, coupled 
with a neutrino interaction probability that scales with 
E, this is equivalent to a 8% uncertainty in flux. For at- 
mospheric neutrinos, the spectrum can be approximated 
by E~ 3 ' 7 for conventional neutrinos and by E~ 2 ' 7 for 
prompt neutrinos from charm and bottom quark decays. 
These lead to uncertainties of 20% and 12%, respectively, 
for the detector sensitivity. We conservatively use the 
larger value for all atmospheric neutrinos. 

Uncertainty in the atmospheric neutrino flux forms an 
additional uncertainty in the atmospheric neutrino back- 



ground to the extraterrestrial cascade search. The ex- 
pected rates of atmospheric neutrinos are based on pre- 
dictions from Ref. [37] and Ref. [7J for the conventional 
and prompt components, respectively. We assign a 20% 
uncertainty in these rates |13) . 

In the relevant energy range, uncertainty in the neu- 
trino interaction cross-sections is about 6%, caused 
largely by uncertainties in the parton densities ED] ■ 

The total uncertainty in the atmospheric (extraterres- 
trial) neutrino detection efficiency was estimated to be 
29% (22%) by adding these contributions in quadrature, 
see Table Hill 



TABLE III: Table of systematic uncertainties for the atmo- 
spheric and extraterrestrial neutrinos. 



Source 


^atm 


u(E- 2 ) 


Ice Properties 


20% 


20% 


DOM Efficiency 


20% 


8% 


Neutrino Cross-Section 


6% 


6% 


Total Uncertainty 


29% 


22% 



To check the energy scale, we studied the detector re- 
sponse to a 337 nm N2 laser, known as the standard 
candle. This laser is on one of the IceCube strings, at a 
depth of 2153 m. It produces light that is shaped like a 
Cherenkov cone pointing downward; the light output is 
calibrated to ±10%. For DOMs that are far enough from 
the laser to avoid saturation (defined as observing less 
than 20,000 photoelectrons), the total observed charge is 
3% lower than our Monte Carlo prediction for the cal- 
ibrated ice model, and 12% lower than the expectation 
for the alternate ice model. This is well within the ex- 
pectations from the ice model and the DOM efficiency 
uncertainty. 
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IV. RESULTS 

After the event selections had been finalized, the anal- 
yses were unblinded. The 10% of the data, which was 
used to develop the selections, was discarded from the 
results. The remaining data correspond to 257 days of 
livetime. 

The effective area A e g, defined as the equivalent area 
with 100% neutrino interaction probability and detec- 
tion efficiency, was obtained by passing simulated sig- 
nal Monte Carlo events through the analyses. Figure [T0| 
shows A e g for v e , , and v T versus neutrino energy. The 
contribution from the Glashow resonance is clearly visi- 
ble in A e g for v e at trigger level and at final selection level 
for the extraterrestrial analysis. The analyses cover com- 
plementary energy ranges, with the atmospheric analysis 
naturally covering lower energies than the extraterrestrial 
analysis. The analyses have similar A g in the region of 
overlap. 

The energy resolution is A(log 10 E v ) ~ 0.26 (0.18) for 
the E~ 2 (atmospheric) v e energy spectrum. The x, y and 
z position resolution is ~ 10 meters. 



A. Atmospheric Neutrino Results 

Before unblinding, the event selection criteria for the 
atmospheric analysis were optimized for a conventional 
atmospheric neutrino flux with a 5 TeV threshold on the 
reconstructed cascade energy. Twelve events remained in 
the 90% of the recorded data after an a-posteriori selec- 
tion criterium on the multivariate selection variable Q* at 
0.9. The expected atmospheric electron and muon neu- 
trino signal is 7.0 ± 2.0 events for these selection criteria. 
The prompt neutrino component is 25% of the expecta- 
tion. Only one standard CORSIKA Monte Carlo event 
passed the final selections. Muon background estimate 
based on this one Monte Carlo event is 9 ± 9. It includes 
an empirical scaling factor to normalize the simulated 
background to the data. The limited background Monte 
Carlo statistics preclude detection of a signal. 

Figure [Tl] shows the distribution of the variable Q* for 
the full data set with -E reC o > 5 TeV and compares it to 
the corresponding distribution for a simulated signal of 
cascades from atmospheric electron and muon neutrinos. 

Figure [12] shows the cumulative number of events as a 
function of the threshold applied on the Q* variable. The 
top panel compares the experimental data to the signal 
expectation from simulations, shown as a band with a 
width determined by the systematic uncertainty in the 



atmospheric neutrino flux. The bottom panel in Fig. 12 
shows the ratio of data events over expected signal events 
as a function of selection threshold. This ratio crosses 
unity for the most discriminating selection criterium that 
can be applied without rejecting all the data, and where 
the remaining few events have the highest probabilities of 
being due to atmospheric neutrino cascades. This implies 
that, in using this analysis, the detector becomes sensi- 
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FIG. 11: Distribution of the multivariate selection parameter 
Q* for data and for a simulated cascade signal from v e and 
with reconstructed energies above 5 TeV in the atmospheric 
analysis. 
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FIG. 12: Experimental events reconstructed with energies 
above 5 TeV in the atmospheric analysis compared to an ex- 
pected cascade signal from atmospheric neutrinos, (top) The 
cumulative number of events above a threshold in the multi- 
variate selection parameter Q* for data and simulated signal, 
(bottom) The ratio of data to expected signal as a function 
of selection threshold. 



tive to atmospheric cascade events due to neutrinos, but 
the exposure is not sufficient for a statistically significant 
detection. 

At a 15 TeV energy threshold and an a-posteriori se- 
lection criterium on the multivariate selection variable 
Q* > 0.9, six events remained in the 90% of the recorded 
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data. The expected background from atmospheric muons 
was estimated to be 5.0 ± 3.8 events. This estimate was 
made with high-energy CORSIKA simulations and in- 
cludes the empirical scaling factor to normalize the simu- 
lated background to the data. The expected atmospheric 
electron and muon neutrino signal is 3.0 ± 0.9 events for 
these selection criteria. The prompt neutrino component 
is 40% of the expectation. The observed number of events 
is compatible with Monte Carlo prediction. 



B. Diffuse Flux Limit 

A total of 14 events passed all selection levels in the ex- 
traterrestrial analysis with an expected total background 
of 8.3 ± 3.6 events. The background was estimated from 
simulations, which indicate sizable contributions from at- 
mospheric muons and from atmospheric neutrinos (see 
Table 

Four events from the high-energy CORSIKA Monte 
Carlo background sample satisfy all event selection cri- 
teria. They all originated from proton induced show- 
ers, with zenith angles of around 70 degrees and ener- 
gies in the range of 0.5 — 3.5 PeV. Figure [15] shows one 
out of these four Monte Carlo muon background events. 
Taking into account the event weights in the simulation 
and the rate normalization, the four simulated events 
that satisfied all selections correspond to 5.4 ± 3.5 muon 
background events in 257 days of livetime. The number 
of background atmospheric neutrinos from conventional 
and prompt sources was estimated to be 2.9 ± 0.9, where 
the central value 2.9 was obtained assuming atmospheric 
neutrino fluxes from Refs. [7J [37] . 

The expected rate of signal neutrino events was 
evaluated for an assumed flux of -E 2< f> m odci = 1-0 x 



10~ b GeV • cm" 



sr for the sum of all flavors and 



using the effective areas, A c g , given in Figure [TU] The 
results are given in Table [TT] for each flavor. In 257 days 
of livetime, a total of 46 ± 10 signal events is expected, 
assuming that this <& m odei receives equal contributions 
from all flavors (v e : : v T flux is 1 : 1 : 1 at the Earth). 
Electron neutrinos contribute about 40%, tau neutrinos 
about 45% and muon neutrinos the remaining 15%. 
The energy distributions are shown in Fig. 
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By using the Feldman-Cousins method [41] as imple- 
mented in the software package POLE++ [42 to ac- 
count for statistical and systematic uncertainties in the 
estimated background and signal counts, we set an up- 
per limit on the number of signal events at 90% confi- 
dence level of /x 90 % = 16.6. This /Kg % is below the ex- 
pected total number of signal events and, hence, $ mo doi 
is rejected at 90% confidence. Conversely, an upper 
limit at 90% confidence can be set at _E 2 $ 



90% C.L. 



< 



3.6 x 10 -7 GeV • cm~ 2 • s _1 • sr _1 on the diffuse flux of 
neutrinos of all flavors assuming that <f> cx E~ 2 and that 
the flux at the detector receives equal contributions from 
all flavors. 

An upper limit on the flux of v e that does not depend 
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FIG. 13: Cascade reconstructed energy distribution after all 
selections in the extraterrestrial analysis. 



on the assumption of equal flux contributions for each 
flavor was derived by assuming that the and v r fluxes 
are zero. This upper limit on the flux of electron neutri- 
nos is £ 2 $ 90 % c.l. < 3.0 x 10~ 7 GeV • cm" 2 • s" 1 • sr" 1 . 

In these limits, the central 90% of v signal events are 
in the energy range from 24 TeV to 6.6 PeV with a mean 
energy of ~ 220 TcV. The limits are shown in Fig. 14 



v 1(T 

</) 
(/) 

CM 

E . 
«10"' 
> 

(5 



>10 



10 ! 



lceCube-22 cascades (all flavor) 257 d 
AMANDA-II cascades (all flavor) 1001 d [13] 
AMANDA-II UHE (all flavor) 457d [43] 
Bniknl 1038 d (all flavor) [44] 



- ■ WB 199BX3/2 [11] 

- - WB Prompl GRB [45] 

- - - Blazars Sleeker 2005 [46] 



9 

/GeV) 



FIG. 14: (Color online) Experimental upper limits on the 
diffuse flux of neutrinos from sources with $ oc E~ 2 energy 
spectrum and theoretical predictions for neutrino fluxes from 
astrophysical neutrino sources. 

Figure [16] shows the experimental event with the high- 
est reconstructed energy that passed all selections in both 
atmospheric and extraterrestrial analyses. 



V. SUMMARY 

In summary, we report the first search for cascades in- 
duced by atmospheric and by diffuse astrophysical neu- 
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FIG. 15: (Color online) The MC muon background event 
which passed all selections in extraterrestrial analyses and 
which has the reconstructed energy of 18 TeV is displayed 
from the side. Different colors of the circles represent differ- 
ent DOMs hit times, with early hits in red and late hits in 
blue. The size of the circles represents the amplitude. 




FIG. 16: (Color online) The data event which passed all selec- 
tions in atmospheric and extraterrestrial analyses and which 
has the highest reconstructed energy of 134 TeV is displayed 
from the side. Different colors of the circles represent different 
DOMs hit times, with early hits in red and late hits in blue. 
The size of the circles represents the amplitude. 



trinos with the IceCube detector. The data, obtained 
in 2007-2008 with a configuration of 22 active strings, 
amount to 257 days of livetime and was searched for 
charged current interactions of v & and v r , and for neutral 
current interactions of neutrinos of all flavors. 

The atmospheric neutrino analysis used neural- 
network based event selections and resulted in a total 
of 12 candidate events with energies above 5 TeV after 
event selections. Within the large uncertainties, the ob- 
served number of events is consistent with the expected 
background. 



The astrophysical neutrino analysis used one and two 
dimensional selection criteria and was optimized for 
higher energies than the atmospheric neutrino analy- 
sis. A total of 14 events with energies above 16 TeV 
remained after event selections, with an expected total 
background contribution of 8.3 ± 3.6 events. We de- 
rive an upper limit at 90% confidence of E 2 $ 90 % CL < 
3.6 x 10~ 7 GeV • cnr 2 • s^ 1 ■ sr _1 on the diffuse flux of 
astrophysical neutrinos with the assumption that the en- 
ergy spectrum <E> ex E~ 2 and that the flavor composition 



of the v P 



v T flux is 1 : 1 : 1 at the Earth. In this 



limit, 90% of the expected signal events have energies 
between 24 TeV and 6.6 PeV. 

This is below the limit that was recently reported from 
final AMANDA data, corresponding to 1001 days of live- 
time [13] ■ Once construction is completed in 2011, Ice- 
Cube will consist of 86 strings covering a volume of 1 km 3 . 
Future IceCube searches will thus benefit from a consid- 
erably larger size and are expected to have significantly 
improved detection sensitivity. 
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